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Abstract 

The long-range interactions of two atoms, of an atom and a dielectric wall, 

of an atom and a perfectly conducting wall, and of an atom between two 

perfectly conducting walls are calculated, including the effects of retardation, 

for Li using dynamic polarizabilities determined from highly correlated, vari- 

ationally determined wave functions. 
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I. INTRODUCTION 



Long-range interactions between two ground state atoms and between a ground state 
atom and a surface are now measured using lasers and cold atoms or atomic beams. Pho- 
toassociation spectroscopy has yielded strict limits on the values of coefficients of dispersion 
forces between two Rb atoms [|T|, two Na atoms 0, and two Li atoms |^ in their ground 
states. Spectroscopy combined with deflection of alkali-metal atomic beams near surfaces |^ 
or reflection of Na atoms from surfaces in atomic fountains |^,^ have made it possible to mea- 
sure the coefficients of atom-surface forces. The experiments are consistent with theoretical 
models for the interaction potentials, but accurate theoretical estimations of the potentials 
remain elusive for atoms other than H and He. The effects of retardation, due to the fi- 
nite speed of light, cause the potentials to become weaker, approaching simple power laws 
for asymptotically large distances. The advent of highly-correlated basis sets for Li using 
multiple nonlinear variational parameters |^ makes it possible to perform well-converged 
calculations of the dynamic electric polarizability functions, thereby enabling, as we will 
show, precise evaluation of long-range interaction potentials, including retardation, for two 
Li atoms, for a Li atom and a dielectric or perfectly conducting wall, and for a Li atom 
between two perfectly conducting walls. 



II. FORMULATION 



A. Atom-atom interactions 



The effect of retardation on the long-range induced dipole interactions of two atoms was 
first investigated by Casimir and Polder p| and the effects on higher induced multipoles by 
Au and Feinberg Jenkins, Salam, and Thirunamachandran and Power and Thiruna- 



machandran [11|. In this paper, the dipolar-dipolar and dipolar-quadrupolar interactions 



are considered, the higher multipolar interactions being negligible. We use the expression for 



the retarded dipolar-quadrupolar interaction of Power and Thirunamachandran W^ which 
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differs from the approximate expression obtained by Au and Feinberg 0. 

The interaction potential between two like atoms including the effects of retardation can 



be written IB. 11 



where 
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and the retardation coefficients are 



where 
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Pii(x) =x^ + 2x-' + 5x^ + 6x + 3 
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and 



where 
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fs{R) = ir~Fr / dujexp{-2aisUjR)ai{iuj)a2{iuj)Pi2{ujaisR), 
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Pi2(x) = + 3x^ + f + 42x3 + 81x2 + 90x + 45, (8) 



and afs = 1/137.035 989 5 is the fine structure constant. We use atomic units throughout. 

The functions ai{iuj) appearing in (0)-(0) are the dynamic electric multipole polarizabil- 
ity functions at imaginary frequency defined by expressions (6)-(9) of ||12|| . 



The retardation coefficients are dimensionless and can be expanded for small R as 

CMR)r^C,-alR'W, (9) 



with 



1 



oo 



W4 = - / dujuj^aiiiuj) (10) 
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and 



CsfsiR) Cs ~ aiR'W,, (11) 



with 
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Wq = — / duj uj'^ai{iuj)a2{iuj). (12) 
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The coefficients W4 and can also be derived from an analysis in perturbation theory 
of the orbit-orbit term arising from the Breit interaction in the Pauli approximation [p!3[ . 



Expanding Power and Thirunamachandran's result for —Csfs{R) for small R, according to 
(|rT|), we find a value of Wq a factor of | times larger than the value of Wq resulting from 
the theory of Au and Feinberg (see Eq. (4.21) of [Q). This resolves the discrepancy of | 
found in between the value of Wll,a;2 for H from Johnson, Epstein, and Meath [[15 



who evaluated terms from the Breit interaction in the Pauli approximation [13|, and the 



value of evaluated for H in |14| using the Au and Feinberg formulation. 



For asymptotically large R, the retardation coefficients have the limits 

23 aj{ 
inafsR C( 
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and 



f(m 531 ai(0)a2(0) „ _ .... 

B 7^ ' R^oo. (14) 
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B. Atom-wall interactions 



Expressions for the interaction potential of an atom and a dielectric wall, including 
the effects of retardation, have been given by Dzyaloshinskii, Lifshitz, and Pitaevskii |]TB 



Parsegian |T^, and Tikochinsky and Spruch [113]. For a wall with a dielectric constant e, the 



potential can be written 



/.oo poo 

VAtD{R,e) = --^ d^ea,{tO dpexp{~2^Rpaf,)H{p,e), (15) 
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where 



H{p,e) = '—^ + {l-2p'y^ (16) 
s + p s + ep 

and 

s = (e-l+p2)i/2. (17) 

We follow the notation of and the subscripts At, D, and M denote, respectively, an 
atom, a dielectric wall, and a perfectly conducting {i.e. metal) wall. 
For asymptotically large distances, 

VAtoiR, e) - V^j,{R, e) = -|ii^0(e), (18) 



where 



and 



K4 = 3ai(0)/(87rafs) = 16.36 ai(0) (19) 



Direct integration of ( |20D yields 



^ e- 1 

where 



1 4-(e + l)ei/2 ^ ^ 



(21) 
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^M=- "'';-;/''' ii+.+2.(.-in (22) 

and 

5(e) = - — ^^[Arcsinh(e^/2^ - Arcsinhfe-^/')], (23) 
(e + 1)^/^ 

in agreement with Dzyaloshinskii et al. [Jl6[. Approximations to V^^(-R, e) will be considered 
in Sec. below. 

The potential for the interaction of an atom and a perfectly conducting wall follows by 



letting e — cxD in ([T5|) giving [I 



where the coefficient is 



Vmm{R) ^ Vmd{R, oo) = (24) 
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Cs = — dua^itcu), (25) 
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and the retardation coefficient is 



1 

/3(^) = B / dx e-^ai(^x/2afei?)[ia;2 + x + 1]. (26) 



Eq. ( p6D approaches for asymptotically large distances the form 



giving 

VAtM{R) - VX^Af (^) = -KJR\ (28) 

The interaction potential for an atom between two parallel, perfectly conducting walls 
has been given by Barton pO| and by Zhou and Spruch It can be expressed as 



VMAtMiz, L) = T2{L) - Ti{z, L), (29) 

where 
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lat^Ls sinh t 

where L is the interwall distance and z is the distance of the atom from the midpoint. For 



small values of L, the potential is [pT 



VMAtM{z,L) -—T{z/L)C3, 



(32) 



where 



JO sinh t 
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and C3 is defined in (^51) . For asymptotically large values of L, the potential is 



7r3ai(0) 



360 



3 - 2cos2(7r^/L) 
8cos4(7rz/L) 



(34) 



III. CALCULATIONS 



The calculations of the wave functions and the polarizability response functions have been 
described previously in, respectively, and [|12[. We briefly summarize the procedures. 



The basis set for the lithium atom is constructed in Hylleraas coordinates 

{<Pt,^, {at, Pt, It) = rf ri' rj^^ ^.3ig-a.ri-/3.r.-^..3| ^ (35) 

where fit denotes a sextuple of integer powers ji, j2, js, ju, j23, and J31, index t labels 
different sets of nonlinear parameters at, (3t and 74. Except for some truncations, all terms 
are included such that 

3i + 32 + is + J12 + J23 + J31 < ^ • (36) 
The wave function is expanded from the multiple basis sets 



^(ri,r2,r3) = AY,Y.^i^t^t<Pt,t,t{(^t, Pult) 



t pit 



X 



(angular function) (spin function) . (37) 



A complete optimization is performed with respect to all the nonlinear parameters. The 
screened hydrogenic wave function is also included explicitly in the basis set. 

The dynamic polarizabilities are evaluated using effective oscillator strengths and tran- 
sition energies obtained from the diagonalization of the Hamiltonian in a basis set of S 
symmetry for the ground state and of P and D symmetry, respectively, for the intermediate 
states corresponding to the dipole and quadrupole polarizabilities. The basis sets were the 
size 919 set from [|12| for the S symmetry and the size 1846 sets from [0 for the P and the 



D symmetries. A detailed discussion of the evaluation of ai{iuj) can be found in |]T2|. The 



static polarizabilities have the values ai(0) = 164.111(2) and 02(0) = 1423.266(5) [1/ 

Values of the coefficient W^^ for two Li atoms have been determined by Margoliash and 
Meath |^ and by Easa and Shukla [^. Using our functions ai{iuj) and 0:2 (^cc'), we de- 



termined the coefficients and using, respectively, (|T0]) and (|T2|), and the results are 
compared with previous results in Table |. We also calculated the coefficients feiR) and 
fsiR) using (H) and (0) at various values of R. The results are given in Table |l| and 
Fig. |l|. The values of the dipole-dipole potential —CQfQ{R)/R^ are in agreement with, but 
are more accurate than, those given in Ref. Jl^, calculated using a model potential method. 



The values of the dipole-quadrupole potential —Csf8{R)/R^ replace those given in Ref. JTi 
which were calculated using the expression of Au and Feinberg as discussed above. The 
dipole-quadrupole potential is usually of secondary importance due to its 1/R^ power law 
behavior. 

Using the polarizability function ai^iu) we evaluated VAtoiR,^) for values of e = 2.123 
and 2.295 corresponding to, respectively, fused silica and BK-7 glass. The values are listed 
in Table ^ and illustrated in Fig. |^ for values of R up to 5000ao. For larger values of R, 
the potential can be obtained from (0). The values of 0(e) from our calculations are listed 
in Table [V| and they are in agreement with the representative values given in Fig. 10 of 
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Ref. m. 



Three approximations V^^^,, V^tD' ^axd ^Ato were obtained by Spruch and 
Tikochinsky by imposing the requirements that the interaction be exact for e ~ oo and for 
e ^ 1, see Eqs. (4.5), (4.9), and (4.12) of Expressing the approximations as ratios to 
the exact potential at very large distances, we have 



T/OO ^ _|_ 37 ' 



23 



V" 



AtD 



F(. 



T/OO 



I 30 1/2 , X' 
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(38) 



(39) 
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and 



T/W 

T/OO 



where 



Fie) 



e + 2 ' 



1 



(40) 



(41) 



We calculated the ratios appearing in (|38|)-(^0D using our values of (j){e) and the results are 
presented in Fig. ^. Our results indicate that the second approximation defined by ( p9D is 
the most accurate, differing by about 6% at most from the exact value of the potential. The 
third approximation was developed for small values of e where it is seen to be somewhat less 
accurate than the second approximation. 

The interaction potential for a Li atom and a perfectly conducting wall was evaluated 
from (P^)-(PB|). The value for the coefficient C3 is in excellent agreement with previous 
determinations, listed in Table ||, particularly with those calculated from the alternative 
expression 



1 

12 







N 
\i=l 







(42) 



which follows from integration of (P^D, where N is the number of electrons. Note that only 
the ground state wave function is required to evaluate (|4^) . For Table || we used expectation 
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values given by King [^] and Yan and Drake [0]. The values of R^V{R) were calculated 
from ( p6D and values are listed in Table |rV|. The present calculations of the potential values 
are in agreement with, but are more accurate than, those given in Ref. . 

The potential VMAthiiz, L) was evaluated using p9| ) for a range of wall separations L 
and distances z of the atom from the midpoint. Values of the energy shift arising from the 
potential for values of L and z that might be realized in an experiment are given in Fig. ^. 

The expressions involving dielectric walls in this paper were obtained under the approx- 
imation e{uj) ^ e(0) = e, where e{uj) is the frequency-dependent dielectric function of the 
wall. This is an excellent approximation for R ~ oo, but at smaller R it could lead to 
significant error if resonances play a role. 
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TABLES 

TABLE L The coefficients W4 and Wq for two Li atoms. Numbers in parentheses represent 
theoretical uncertainty due to the finite basis set size. 

W4 We Reference 

3.214(2) 219.9(2) Present 



2.9312 Easa and Shukla [|3| 



3.233 Margohash and Meath |22 
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TABLE II. The coefficient C3 for the Li atom-wah interaction calculated in the present work 
compared to values calculated from (^), (direct), using matrix elements given by various authors, 
or from pseudo oscillator strength distribution data of dimension M tabulated by various authors 
(osc. str.). 

C3 Method Reference 

1.518(2) osc. str. Present 

1.518 000 51(3) direct Yan and Drake 

1.518 000 direct King H 



1.49 osc. str. (M = 11) Stacey and Dalgarno [26| 



1.52 osc. str. (M = 10) Margohash and Meath |22 
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TABLE III. The dimensionless retardation coefficients fe{R) and fsiR) for the atom-atom 
interaction. The dispersion coefficients Cq and Cg from [l^ are also given. 





1393.39(16) 


83425.8(4.2) 


R 


feiR) 


fsiR) 


15 


1.0000 


1.0000 


20 


1.0000 


0.9999 


25 


0.9999 


0.9999 


30 


0.9999 


0.9999 


50 


0.9997 


0.9997 


70 


0.9995 


0.9994 


100 


0.9991 


0.9988 


150 


0.9980 


0.9974 


200 


0.9966 


0.9955 


250 


0.9950 


0.9933 


300 


0.9931 


0.9907 


500 


0.9833 


0.9775 


700 


0.9708 


0.9608 


1000 


0.9489 


0.9319 


1500 


0.9076 


0.8791 


2000 


0.8641 


0.8256 


2500 


0.8208 


0.7743 


3000 


0.7789 


0.7263 


5000 


0.6341 


0.5709 


7000 


0.5253 


0.4627 


10000 


0.4113 


0.3555 


15000 


0.2970 


0.2528 
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20000 
25000 
30000 
50000 
70000 
100000 



0.2304 
0.1875 
0.1578 
0.0961 
0.0689 
0.0484 



0.1947 
0.1579 
0.1326 
0.0805 
0.0577 
0.0405 
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TABLE IV. Values of —R^VAtD{R, e), where VAtoiR, e) is the atom-wall potential, for values 



of e corresponding to fused silica and BK-7 glass 


in, respectively, cols. 2 and 3, and in 


col. 4 values 


of -R^VAtM{R) for 


a perfectly conducting wall. 








Fused silica 


BK-7 glass 


perfect 


R 


e = 2.123 


e = 2.295 


e = oo 


10 


0.5360 


0.5859 


1.5007 


15 


0.5323 


0.5819 


1.4937 


20 


0.5289 


0.5782 


1.4871 


25 


0.5259 


0.5749 


1.4810 


30 


0.5230 


0.5717 


1.4753 


50 


0.5130 


0.5608 


1.4551 


70 


0.5045 


0.5515 


1.4380 


100 


0.4933 


0.5392 


1.4157 


150 


0.4772 


0.5215 


1.3836 


200 


0.4629 


0.5060 


1.3551 


250 


0.4500 


0.4919 


1.3289 


300 


0.4381 


0.4788 


1.3042 


500 


0.3974 


0.4344 


1.2160 


700 


0.3644 


0.3983 


1.1393 


1000 


0.3244 


0.3546 


1.0398 


1500 


0.2741 


0.2996 


0.9049 


2000 


0.2368 


0.2589 


0.7981 


2500 


0.2081 


0.2276 


0.7118 


3000 


0.1853 


0.2026 


0.6409 


5000 


0.1276 


0.1395 


0.4526 
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TABLE V. The dimensionless function ^(e). 



e 


He) 


1 


23 
30 


1.5 


0.761364 


2 


0.760757 


2.123 


0.760970 


2.295 


0.761425 


4 


0.770171 


7 


0.787334 


9 


0.797062 


13 


0.812791 


16 


0.822186 


20 


0.832501 


50 


0.874337 


100 


0.902534 


500 


0.950261 


1000 


0.963647 


5000 


0.982986 


10000 


0.987836 


50000 


0.994478 


1 X 10^^ 


0.999996 
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FIGURES 
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FIG. 2. Values of R'^VAtoiR, e), for values of e eorresponding to fused silica (e = 2.123), BK-7 
jlass (e = 2.295), and for a perfectly conducting wall (e = oo). 
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FIG. 3. Ratio of several approximations given by Spruch and Tikochinsky to the the exact 
atom-wah potential calculated in the present work. The symbols A, B, and C represent, respec- 



tively, values from Eqs. (gSD, (|39D, and (|40D. 
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FIG. 4. The energy shift arising from the wall-atom-wall potential for various values of the wall 
separation L and the distance of the atom from the midpoint z. Only the values for z > are 
shown as the energy shift is symmetric about the z = plane. 
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